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Abstract 

The Bar khausen jumps or avalanches in magnetic domain-walls motion between succesive pinned 
configurations, due the competition among magnetic external driving force and substrum quenched 
disorder, appear in bulk materials and thin films. We introduce a model based in rules for the 
domain wall evolution of ferromagnetic media with exchange or short-range interactions, that 
include disorder and driving force effects. We simulate in 2-dimensions with Monte Carlo dynamics, 
calculate numerically distributions of sizes and durations of the jumps and find power-law critical 
behavior. The avalanche-size exponent is in excellent agreement with experimental results for thin 
films and is close to predictions of the other models, such as like random-field and random-bond 
disorder, or functional renormalization group. The model allows us to review current issues in 
the study of avalanches motion of the magnetic domain walls in thin films with ferromagnetic 
interactions and opens a new approach to describe these materials with dipolar or long-range 
interactions^ 
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I. INTRODUCTION 



Depending on shape, composition, and energy contributions, among others, ferromag- 
netic materials can form magnetic domains. These domains are separated by domain walls 
(DW) and their magnetization is slightly inhomogeneous, due to internal disorder. The 
disorder includes vacancies, defects, impurities and dislocations, fixed inside the material, 
and limits the movement of DW when the material is placed in an external magnetic field. 
The quenched disorder, known as Barkhausen noise by the name of the scientist who first 
observed, jl] is an early example of the origen of complexity in materials science and has been 
intensively studied to date.|2|,|3| In many cases, the wall exhibit a jerky motion between states 
of quiescence, known as Barkhausen jumps or avalanches. The competition between restor- 
ing forces and quenched disorder results in transitory mult i-metast able states. Locally, the 
wall leaves a metastable state basically by two mechanisms: by fluctuations spontaneously 
or, under certain conditions, driven by magnetic external force. Examples of avalanches 
with quenched sustrate disorder are the dislocation and crack propagation, (4] stick-slip mo- 



tion of tectonic plates 



5| and those with origin from depinning of fluid contact-interface. {(J 
Avalanches have also been studied in models without quenched substrate disorder, such 
as in sandpile models and in granular matter. [7j An important characteristics of avalanche 
motion is its scale invariance, self-organized criticality, a distribution P(S) ~ S~ T of the 
avalanche- sizes S, and a distribution P(T) ~ T~ a of the avalanches-durations T, being r 
and a the size and duration critical exponents, respectivelly. 



The structure and motion of the DW in bulk materials with defects has been studied long 



time. 



8j In such systems the more accepted structure is that of Bloch, where the magnetiza- 



tion rotates through of wall, outside of the plane of the magnetic moments. The statistical 
properties of the DW avalanches are understood today by means of depinning transition 
models, 0, [hJ] that explain very well many experimental results. Models and experiments 
have been grouped into two universality classes depending on the values of critical expo- 



nents, which also verified scaling equalities fully accepted, llj One class includes materials 



with dipolar or long-range interactions, e.g. polycrystaline and crystallized amorphous ma- 
terials, where the exponents take values close to T£ = 1.5 and ae = 2.0. [12] The other 
class includes materials with predominantly exchange or short-range interactions, e.g. soft- 
magnetic materials with high anisotropy, where the exponents take values close to r s = 1.27 



and a s = 1.50 .0, Q, Q When dipolar interactions are neglegect front propagation occurs, 
similarly e.g. fluid imbibition in porous media, otherwise there is a mean field like behavior. 



In thin films the wall motion is dominated by depinning too, but additionally the DW 
change their structure depending on the thickness of the film. Neel was the first to note that 
the theory of the Bloch walls is not valid for thin films if the thickness is comparable to the 
width of the wall. jl5| As the thickness decreases, the wall passes through the stages of Bloch 
wall, asymmetric walls, crosstie wall, and Neel wall where atomic moments remain in plane 
while the rotation occurs. [3] It is not known yet if the wall structure affects movement. 
Some experimental results show that no changes expected in the statistical results when 
varying the thickness of the thin film. 17] However, theoretically, the manner of interaction 
of the wall with the defects could determine whether the magnetization can rotate locally 



overcome the defect or, in other words, if the wall can pass across the defect. [18] The thin 
film exponents should be universal in the sense of describing a large number of materials and 
systems and must be independent of microscopic details. Magnetic thin films (bellow 200 
nm) have exponents belonging to a new universality classes different from the 3-dimensional 



magnets. 



19H21j The usual theoretical approach to these thin films are 2-dimensional mod- 



els. Recently, experimental results show that the critical exponent of the size distribution 
decreases with increasing temperature, a fact attributed to a transition between the uni- 



versality classes of long- to short-range interactions. 



2l| However, theoretically, this is an 



argument that can be sustained only for 3-dimensional systems. 22] 



The paper is organized as follows. In the next Section we present the simulation model 
and the statistical properties of the DW avalanches. We introduce a model based on local 
evolution rules of the wall, with or without external magnetic field, taking into account 
the defects found the wall in your progress. Also we present the definitions of avalanche 
size and duration used here to calculate the distributions. In Section III we provide an 
overview of our numerical results. We present the basic statistical properties of the distribu- 
tions of avalanche-size and -duration. Also we analyze the results compared with available 
experimental and theoretical. Section IV finishes the paper with conclusions. 
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II. THE MODEL 



In a magnetizable material we consider a volume element consisting of two domains 
separated by a magnetic DW. We assume both domains composed of opposite magnetic 
dipoles and point defects arranged in the nodes or sites on a regular array. We call active 
(or inactive) to sites containing dipoles (or defects). We suppose that the DW consists 
only of active sites and can move through the inactive ones. In the model we disregard the 
possible mechanisms that explain the structure of the DW. Just assume that at any moment 
the DW is a monolayer dipole with magnetic dipole moment perpendicular to the domains, 
and that the magnetization is in the "easy direction" . 

To model a thin film, we consider dipoles and defects located at the nodes of a two- 
dimensional square lattice of edge L with periodic boundary condition. In this lattice, we 
assign a random pinning function k) uniformly distributed in the interval [0, 1] to every 
node (i,k) of the square lattice. For a given parameter p G [0,1], if k) < p the node 
(i, k) is active, otherwise is inactive. To characterize the disorder in the lattice we use the 



activity function F(i,k) = 9(p — £(L fc)).[23| where 9{x) is the unit step function defined 
as 9{x) = 1 for x > and 9{x) = for x < 0. A dipole located at the nodes of the DW 
has perpendicular component to the easy direction. If (j,rij) is a point on the DW in the 
discrete lattice, with j = 1, . . . , L and rij = 1,2,..., an active site (or inactive) next to the 
wall is described by F(j, rij ± 1) = 1 (or 0). 

Our model is based on local evolution rules for all points of the DW. These rules include 
the short-range ferromagnetic interactions, taking into account the local balance of magnetic 
moment on both sides of a point on the wall. When there is local unbalance in the opposite 
direction of movement, the DW at this point spontaneously reaches equilibrium. Otherwise, 
when there is balance, the DW can move driven by the applied field. Hence we assume that a 
point on the DW can move spontaneously with probability c if the magnetic dipole moment 
in the environment is unbalanced on both sides of the wall. Else, if there is an approximate 
balance on both sides, with probability 1— c the DW point on the wall is driven by action of 
the field that attempts to remove it from the metastable state. We start the simulation with 
a flat DW. Chosen column j, the point (j,rij) from the wall evolves according the following 
rules. 

• With probability c : if F(j, rij + 1) = 1 and (a) \{rij+i + rij-x) — rij > 1 the wall point 
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moves one unit [i.e. rij — > rij + 1], or (b) otherwise the wall point is pinned. 
• Else, with probability 1 — c: 

— If F(j, rij + 1) = 1 and rij +1 = rij-\ > rij the wall point moves one unit up from 
its neighbors [i.e. rij — > + 1], otherwise the wall point moves (or not) up to 
the maximum between the site or its neighbors [i.e. rij — > max(n J+1 , Wj-i, n>j)]', 

— If F(j, rij + 1) = and (a) = rij-i > + 1 the wall point moves one unit up 
from its neighbors [i.e. rij — > rij+i + 1], or (b) n,- + i 7^ n^-i and max(n 3+ i, %-i) > 
nj + 2 the wall point moves up to the maximum between its neighbors [i.e. rij — > 
max(nj + i,7ij_i)], or (c) otherwise the wall point is pinned. 

As the DW is made up of active sites, a rule is frustrated when a wall point trying to reach 
inactive sites. 

If n(t) = (ni(t), . . . ,n L (t)) is the position of the DW at time t, it is in metastable state 
during a time interval At > if n(t) = n(t + St) forall 5t < At. Otherwise, the DW moves 
to the next metastable state after a time interval T, making a jump S = X]j=i[ n j(^ + T) — 
rij(t)]. Note that the jump is proportional to the average displacement of the DW between 
metastable states. To calculate the distributions of size S and duration T we extract the 
sequences of avalanches, over many realizations, between first- and last-jump, where the DW 
is finally pinned. 

III. RESULTS 

The distributions of avalanche-size P(S) and avalanche- duration P(T) have been ob- 
tained for different values of c and p and different values os L between 64 and 4096. Powers 
law-type behavior was observed on both distributions. 

In Figure 1 we show P(S) for different densities of active sites p close to criticality 
(Pc — 0.865 ± 0.005) and c fix close to zero. The critical exponent of avalanche size is 
t = 1.29 ± 0.02. It is possible to note that as increases the value of p diminishes the value 
of r. This behavior was observed for different values of L. The inset (b) shows that as 
decreases the value of p the system moves away from the criticality, reducing the region that 
presents a behavior power law type. 
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FIG. 1. (color online) Log-Log plot of normalized avalanche size distribution P(S) with L = 1024, 
c = 10 -3 and different values of p [0.860 (circle), 0.865 (box), 0.870 (diamond)]. Inset: (a) 
Log-Log plot of the data as is were obtained (p = 0.860) and the points considered for a better 
visualization logarithmically equispaced. (b) Log-Log plot of normalized avalanche size distribution 
with different values of p [0.870 (diamond), 0.845 (triangle up)]. 



In Figure 2 we show P(S) for different values of c and density of active sites fix p = 0.860. 
In all the cases, also for different sizes L and density p in the range 0.835 — 0.880, it could 
observe a diminution of r with the increase of c. For c = the system moves away from 
the power law type behavior dramatically as can observe on the inset (a), hence c = can 
be take as another critical point. From inset (b), it is possible to notice that the size of 
the system does not affect considerably the behavior except for very high values of S where 



finite size effect can be observed for the case of lattice with small edge L. 

Our results for r are in excellent agreement with experimental results reported for fer- 
romagnetic thin films. Because the determination of avalanche sizes in thin films is too 
difficult, the amount of experimental data is limited. Recently some data using magneto 
optics measurements was reported. For Co thin films, Kim et al. reported power law behav- 
ior of the avalanche- size distribution, with critical exponent for various film thickness, with 



values between r = 1.34 ± 0.07 (for 5 nm) and r = 1.30 ± 0.05 (for 50 nm) 



20J In a recent 



work Shin et al. reported similar experimental results on Co films, with r ~ 1.33 in the 



thickness range 5-50 nm. 24J For MnAs films Ryu et al. founded that the scaling exponent 
r varies continuosly from 1.32 ± 0.06 to 1.04 ± 0.05 as increase the temperature from 20 °C 
to 35 C.[2l| Magni et al. performed multiscale measurements of avalanche distributions 
varying the observation window size on thin films of permalloy using a hight resolution Kerr 
magnetometer. 25j They determined r ~ 1.18 measuring at several window widths of an 170 



nm film. More recently, Papanikolaou et al. report measures of Barkhausen noise time series 



in 1 /im-thin films of permaloy using techniques inductive in an open magnetic circuit. [26] 
According to our knowledge for ferromagnetic thin films, they first established the power 
law behavior for the distribution of avalanche duration with values for the exponent a in 
the mean-field universality class for 3-dimensional magnets. 

Rosso et al. recently calculated numerically the jumps distribution between succes- 
sively pinned configurations in a 2-dimensional interface spring model of small intensity 



in a random-field landscape. [27J They study mostly random-field disorder (RF), but we also 



check that the results are the same for random-bond disorder (RB) and compare with pre- 



dictions of the Functional Renormalization Group (FRG) field theory. 28j For RF disorder 
they determine r = 1.30 ± 0.02 near to FRG prediction r = 1.2735 ± 0.0005. Although that 
our value for r is close, there are no other numerical or theoretical results for the exponent 
duration a for another comparison. 

Figure 3 shows P{T) for c = 10~ 3 , p = 0.865 and different values of L. It is possible to 
note that on short duration jumps regime all the curves show the typical power law, with 
a = 1.55 ± 0.05. This behavior is limited by the lattice size L. After certain jump duration 
T* for each L the distribution presented a different behavior. This crossover from one regime 
to another could be attribute to finite size effect, wich is consistent with the fact that short 
(or large) duration jumps are possible in small (or giant) lattices. For T > T* the duration 
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FIG. 2. (color online) Log-Log plot of normalized avalanche size distribution with L = 1024, 
p = 0.860 and different values of c [10 -3 (circle), 10~ 2 (box), 10 _1 (diamond)]. Inset: (a) Log-Log 
plot of normalized avalanche size distribution with L = 1024, p = 0.860 and different values of 
c [10 -3 (circle), (triangle up)] (b) Log-Log plot of normalized avalanche size distribution with 
c = 10 -3 , p = 0.860 and different values of L [64 (circle), 4096 (triangle up)]. 



distribution P(T) oc e _T//T ° where was obtained from the simulation that T oc L A with A 
between 1.0 and 1.15. When L — >■ oo is interesting to note that P(T) — > and T* — > oo, 
evidence of finite size effects. While the scaling regime remains essentially unchanged a peak 
occurs around the cutoff. The inset (b) shows in more detail the exponential regime where 
the dashed lines represents the fitting for each curve. 
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FIG. 3. (color online) Log-Log plot of normalized avalanche duration distribution with c = 10~ 3 
and p = 0.865 and different values of lattice edge L between 128 and 1024. Inset: (a) Log-Log 
plot of the data as is were obtained (L = 128) and the points considered for a better visualization 
logarithmically equispaced. (b) Semi-Log plot, detail of the central plot showing the effects of finite 
size, with an exponential fit in dotted line. 

IV. CONCLUSIONS 

The main achievement of our work has been to introduce an avalanches model that 
describes the motion in disordered media of DW for magnetic thin film in which only short- 
range interactions are considered. The simulation results show the power law like critical 
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behavior following the distributions of size and duration of Barkhausen jumps. From the 
results of the simulation, around the criticality, we observe that r decreases when p or c 
increases. At criticality the results for the exponent r of the size distribution is in very good 
agreement with experimental and theoretical work scarce. On the contrary, as the expo- 
nent a of the duration distribution is not been measured experimentally for 2-dimensional 
magnets, our results can not be checked. Comparing with the experimental values for bulk 
magnetic materials where short-range interactions dominate, our values for r are very close 
while those for a differ considerably. Our understanding of two-dimensional modeling for 
ferromagnetic thin-film remains incomplete. Other exponents should be calculated in order 
to establish the scaling relations, e.g. exponent of growth and dynamic exponent, among 
others. In our model the avalanche size distribution does not depend significantly on the lat- 
tice size, since the power-law behavior extends over two logarithmic orders. By contrast, our 
model shows that the avalanches duration distribution includes finite size effects, which tend 
to disappear by increasing the lattice size according to an exponential cutoff. An outstanding 
question in thin films is the theoretical and experimental characterization of materials with 
predominantly long-range interactions. 29| Finally, models based on evolution rules could be 
ussefull to take into account long-range interactions. 



[1] H. Barkhausen, Z. Phys. 20, 401 (1919) 

[2] J. S. Urbach, R. C. Madison, and J. T. Markert, Phys. Rev. Lett. 75, 276 (1995) 

[3] K. Dahmen, S. Kartha, J. A. Krumhansl, B. W. Roberts, J. P. Sethna, and J. D. Shore, J. 

Appl. Phys. 75, 5946 (1994)0. Perkovic, K. Dahmen, and J. P. Sethna, Phys. Rev. Lett. 75, 

4528 (1995)K. A. Dahmen, J. P. Sethna, and O. Perkovic, IEEE T. Magn. 36, 3150 (2000)J. H. 

Carpenter, K. A. Dahmen, J. P. Sethna, G. Friedman, S. Loverde, and A. Vanderveld, J. Appl. 

Phys. 89, 6799 (2001)K. A. Dahmen, J. P. Sethna, M. C. Kuntz, and O. Perkovic, J. Magn. 

Magn. Mater. 226-230, 1287 (2001)A. P. Mehta, A. C. Mills, K. A. Dahmen, and J. P. Sethna, 

Phys. Rev. E 65, 046139 (2002) 
[4] P. Moretti, M. C. Miguel, M. Zaiser, and S. Zapperi, Phys. Rev. B 69, 214103 (2004)L. Ponson, 

D. Bonamy, and E. Bouchaud, Phys. Rev. Lett. 96, 035506 (2006)L. Ponson, Ann. Phys. 32, 

1 (2007)D. Bonamy, S. Santucci, and L. Ponson, Phys. Rev. Lett. 101, 045501 (2008) 

10 



[5] D. S. Fisher, K. Dahmen, S. Ramanathan, and Y. Ben-Zion, Phys. Rev. Lett. 78, 4885 
(1997)D. S. Fisher, Phys. Rep. 31, 113 (1998)J. M. Schwarz and D. S. Fisher, Phys. Rev. 
E 67, 021603 (2003)E. A. Jagla and A. B. Kolton, J. Geophys. Res. B 115, 05312 (2010) 

[6] A. Prevost, E. Rolley, and C. Guthmann, Phys. Rev. B 65, 064517 (2002)Phys. Rev. Lett. 83, 
348 (1999)S. Moulinet, C. Guthmann, and E. Rolley, Eur. Phys. J. E 8, 437 (2002)S. Moulinet, 
A. Rosso, W. Krauth, and E. Rolley, Phys. Rev. E 69, 035103(R) (2004) 

[7] P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987)C. Tang and P. Bak, 
J. Stat. Phys. 51, 797 (1988)D. Dhar and R. Ramaswamy, Phys. Rev. Lett. 63, 1659 
(1989)S. Banerjee, S. B. Santra, and I. Bose, Z. Phys. B 96, 571 (1995)D. Dhar, Physica 
A 369, 29 (2006) 

[8] G. Bertotti, Hysteresis in Magnetism (Academic Press, San Diego CA, 1998) 
[9] G. Durin and S. Zapperi, in The Science of Hysteresis, Vol. II, edited by G. Bertotti and 
I. Mayergoyz (Elsevier, Amsterdam, 2005) Chap. 3, pp. 181-267 
[10] J. P. Sethna, K. A. Dahmen, and O. Perkovic, in The Science of Hysteresis, Vol. II, edited by 

G. Bertotti and I. Mayergoyz (Elsevier, Amsterdam, 2005) Chap. 2, pp. 107-179 
[11] F. Colaiori, Adv. in Phys. 57, 287 (2008) 
[12] G. Durin and S. Zapperi, Phys. Rev. Lett. 84, 4705 (2000) 

[13] M. Bahiana, B. Koiller, S. L. A. de Queiroz, J. C. Denardin, and R. L. Sommer, Phys. Rev. 

E 59, 3884 (1999) 
[14] G. Durin and S. Zapperi, J. Appl. Phys. 85, 5196 (1999) 

[15] L. Neel, C. R. Acad. Sci. Paris 241, 533 (1955)E. E. Huber Jr., D. O. Smith, and J. B. 

Goodenough, J. Appl. Phys. 29, 294 (1958) 
[16] A. Hubert and R. Schafer, Magnetic Domains (Springer, Berlin, 1998) 

[17] L. Santi, F. Bohn, A. D. C. Viegas, G. Durin, A. Magni, R. Bonin, S. Zapperi, and R. L. 

Sommer, Physica B 384, 144 (2006) 
[18] Barkhausen avalanches in Co thin films have been observed directly by D.-H. Kim, S.-B. Choe, 

and S.-C. Shin, J. Appl. Phys. 93, 6565 (2003) 
[19] E. Puppin, Phys. Rev. Lett. 84, 5415 (2000) 

[20] D.-H. Kim, S.-B. Choe, and S.-C. Shin, Phys. Rev. Lett. 90, 087203 (2003) 
[21] K.-S. Ryu, H. Akinaga, and S.-C. Shin, Nat. Phys. 3, 547 (2007) 
[22] J. P. Sethna, Nat. Phys. 3, 518 (2007) 

11 



[23] L. A. Braunstein, R. C. Buceta, C. D. Archubi, and G. Costanza, Phys. Rev. E 62, 3920 

(2000)G. Costanza, Physica A 330, 421 (2003) 
[24] S.-C. Shin, K.-S. Ryu, D.-H. Kim, and H. Akinaga, J. Appl. Phys. 103, 07D907 (2008) 
[25] A. Magni, G. Durin, S. Zapped, and J. P. Sethna, J. Stat. Mech. 01, P01020 (2009) 
[26] S. Papanikolaou, F. Bohn, R. L. Sommer, G. Durin, S. Zapperi, and J. P. Sethna, Nat. Phys. 

7, 316 (2011) 

[27] A. Rosso, P. Le Doussal, and K. J. Wiese, Phys. Rev. B 80, 144204 (2009) 

[28] O. Narayan and D. S. Fisher, Phys. Rev. B 48, 7030 (1993)H. Leschhorn, T. Nattermann, 

S. Stepanow, and L. H. Tang, Ann. Phys. 509, 1 (1997) 
[29] A. Mughal, L. Laurson, G. Durin, and S. Zapperi, IEEE T. Magn. 46, 228 (2010) 



12 



